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Recently it was shown that vacuum in a background of strong enough magnetic field becomes 
an electromagnetic superconductor due to interplay between strong and electromagnetic forces. 
The superconducting ground state of the vacuum is associated with a spontaneous emergence of 
quark- ant iquark condensates which carry quantum numbers of charged p mesons. The p-meson 
condensate is an inhomogeneous structure made of the so-called p vortices, which are parallel to the 
magnetic field axis. The condensation of the charged p mesons induces a (much weaker) superfluid- 
like condensate with quantum numbers of the neutral p'^^^ mesons. In this paper we show that the 
vortices in the superconducting condensate organize themselves in an equilateral triangular lattice 
similarly to an ordinary type-II superconductor. We show that each of these superconductor vortices 
is accompanied by three superfiuid vortices and three superfiuid antivortices made of the neutral 
p meson condensate. The superconductor vortex overlaps with one of the superfiuid vortices. The 
superposition of the superconducting and superfiuid vortex lattices has a honeycomb pattern. 

PACS numbers: 12.38.-t, 13.40.-f, 74.90.+n 



I. INTRODUCTION 

The behavior of quantum field theories in an extremely 
strong external magnetic field of hadronic scale has at- 
tracted growing interest of the scientific community. The 
chiral magnetic effect [1 , which may take place in hot 
quark matter created in heavy- ion collisions [2 , provides 
a particularly interesting and experimentally observable 
example: chirally-imbalanced matter generates electric 
current along the axis of the magnetic field [3 . The inter- 
est in the strong magnetic fields is far from being purely 
academic since extremely high magnetic fields may be 
generated in noncentral collisions of heavy ions [TJ |4], 
and they may have existed in the early moments of our 
Universe [5^. 

Not only (dense) matter, but also a quantum vacuum 
may exhibit quite unusual properties in a sufficiently 
strong magnetic field background. A well-known related 
example is the magnetic catalysis [BHH], which implies, 
in particular, a steady enhancement of the chiral sym- 
metry breaking in the vacuum of Quantum Chromody- 
namics (QCD) as the external magnetic field strengthens. 
Consequently, a strong background magnetic field affects 
drastically the phase structure of QCD vacuum [9HTT] 
and QCD matter [H [13] 

More recently it was shown that a sufficiently strong 
magnetic field of hadronic scale may cause the vacuum to 
behave as an inhomogeneous and anisotropic electromag- 
netic superconductor [Ml [15] . The superconductivity of, 
basically, empty space, is caused by a spontaneous cre- 
ation of a (charged) p-meson condensate if the strength 
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of the magnetic field exceeds the critical value 

Be ^ 10^^ Tesla or eB^ - 0.6 GeV^ . (1) 

The vacuum superconductivity is accompanied by a 
superfiuid- like condensation of the neutral p mesons J4] . 

The charged p mesons - or, better to say, quark- 
antiquark condensates with the quantum numbers of the 
p mesons - play a central role in the superconducting 
mechanism. The p mesons are vector particles with 
anomalously high magnetic moment corresponding to the 
gyromagnetic ratio g = 2. One can argue that the 
anomalous magnetic moment provides a large negative 
contribution to the squared energy of the p mesons in a 
background of strong magnetic fields. As the external 
magnetic field exceeds the critical value of the magnetic 
field ([1]), the energy becomes purely imaginary indicat- 
ing a condensation of the p mesons. The emergence of 
the electrically charged condensate implies, almost in- 
evitably, electromagnetic superconductivity of the new 
vacuum ground state. [14 

One can also argue that in the background of the 
strong magnetic field the charged vector mesons play a 
role of the Cooper pairs [16]: the strong magnetic field 
makes the motion of the quarks essentially one dimen- 
sional because the electrically charged quarks may move 
only along the magnetic lines. In one spatial dimension 
a weak attraction between a quark (for example, "up" 
quark) and an antiquark (say, "down" antiquark) me- 
diated by a (virtual) gluon inevitably leads to creation 
of a bound state, electrically charged vector meson (in 
our example, it is p+ = ud meson). The emergence of 
the bound states leads to lowering of the vacuum energy 
and, again, to condensation of the charged p mesons. 

The superconductivity of the vacuum in a strong 
magnetic field was first found in an effective bosonic 
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model which describes the electrodynamics of the p 
mesons p^. Later, the superconductivity effect was con- 
firmed in the Nambu-Jona-Lasinio model [15]. Signa- 
tures of this counterintuitive effect were also found in 
holographic approaches [17,18^ and in numerical simula- 
tions of quenched lattice QCD [19]. 

Due to the anisotropic nature of the superconductiv- 
ity (the vacuum superconducts only along the axis of the 
magnetic field) the Meissner effect is absent so that the p- 
meson condensate does not screen the external magnetic 
field [m [16]. Moreover, due to the anisotropic super- 
conductivity the vacuum becomes a (hyperbolic) meta- 
material which, electromagnet ically, behaves similarly to 
diffractionless "perfect lenses" [20] . 

This paper is devoted to a detailed study of the vortex 
contents of the superconducting ground state of the vac- 
uum in a background of a strong magnetic field. In the 
superconducting phase the p-meson condensate forms an 
inhomogeneous periodic (lattice) structure made of the 
so-called p vortices which are parallel to the magnetic 
field axis [TT . The p vortex - which is a string-like topo- 
logical defect in the p-meson condensate - is a close ana- 
logue of the Abrikosov vortex in an ordinary supercon- 
ductor. In Section [IT] we briefly outline a few basic prop- 
erties of the Abrikosov vortex lattice in a mixed state of a 
type-II superconductor in the Ginzburg-Landau (GL) ap- 
proach to superconductivity. Following this example, in 
Section [ITT] we discuss various properties of the p-vortex 
lattice such as the geometrical structure, inhomogeneous 
superconductivity and the structure of inhomogeneities - 
including the vortex content - of the neutral (superfluid) 
p-meson condensate. The last Section is devoted to our 
conclusions. 



II. GINZBURG-LANDAU MODEL 

Consider the Abelian Higgs model which is a relativis- 
tic version of the GL model 

£(0, A) = - Jf^.F^-- + - V{<P) , (2) 

V{<P) = -m'\<p''\ + -^\<P\\ 

where = — ieA^ is the covariant derivative, 
is the electromagnetic field with the strength F^^^ = 
d^Ai, — djyA^^ and cj) = e*'^ is the electrically charged 
scalar field with the mass parameter m and the self- 
interaction coupling A. The field (j) - which carries a 
unit (e) electric charge - plays a role of a field of Cooper 
pairs (without loss of generality, we consider the singly- 
charged bosons (j) instead of the usual doubly- charged 
bosons corresponding to the Cooper pairs of electrons). 

In the condensed state, > 0, the ground state of the 
model ([2| is characterized by a condensate {<p) = 0o with 
1 00 1 = Y^2m^/A, while the photon and the scalar 
excitation S(j) = (j) — {(j)) acquire the following masses. 



respectively: 

/- 2em 
= V 2m , rriA = • (3) 
V A 

We introduce an external magnetic field parallel to the 
z axis, which is described by a "symmetric" gauge po- 
tential = A^^S 

Ai=-^X2, ^2 = ^X1, A3=Ao=0, (4) 

and ignore quantum fluctuations of the gauge fields and 
scalar fields, thus treating the problem at the classical 
level. We always assume that ^ext > (and e = |e|) for 
the sake of simplicity. 

We consider a type-II superconductor which corre- 
sponds to the following region of the parameters: 

X>e\ (5) 

If the external magnetic field exceeds a certain (first) 
critical field, 5ext > ^ then the Abrikosov vortices are 
formed in the superconducting material. An elementary 
Abrikosov vortex carries the magnetic flux 

J d?xB{x) = ^, (6) 

where the integration goes over the plane which is per- 
pendicular to the vortex axis. 

In the type-II regime the vortices repel each other 
and they form a regular periodic structure called "the 
Abrikosov lattice" . If the strength of the magnetic field 
exceeds the second critical field, ^ext > with 

Be = Bc^ = -— , (7) 

then the superconductivity gets destroyed completely. 

We work in the regime quite close to the critical 
field 0: 

^ext < \Bc — 5ext| <^ • (8) 

The ground state is independent of the time t and z coor- 
dinates, so that the equations of motion of the system ([2| 
can be written in the complexified form: 

S(5 + e|0p) -2e0^D0 = O, (9) 
(DI) + - eB)(t) - A|0p0 = , (10) 

with z = z — Xi—ix2, d = 81^182-, d = 81 — 182, 

A = ^1+^2, A = Ai-iA2, B = Fi2 = -i{dA-8A)/2, 
and the covariant derivatives are: 

D = 5 + ^^ext^ , I) = a + ^^ext^ , (11) 

As the magnetic field increases towards B^ the con- 
densate diminishes and at the critical value of the field 
the condensate vanishes completely, 0(5ext = ^c) = 0- 
Thus in the regime (|8|, the classical equations of motion 
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^ and ([To]) can be linearized. In the leading order we 
get the following equation for the condensate (j): 



The solution (j) 
magnetic field, 



D^ = 0. (12) 
{z) of this equation determines the 



B{z) = B, 



(13) 



which is a certain function of the transversal coordinates 
xi and X2- The brackets (...) indicate a mean value in 
the transversal (xi,X2) plane: 



{O) 



Area 



^ — / dxi I 
ea^ J J 



dx20{xi,X2), (14) 



where Area^ is the area of the transversal plane. An 
additive coordinate-independent term in the so- 

lution (13) is a normalization term which imposes the 



conservation of the magnetic fiux coming through the 
transversal plane, 



{B) = Be: 



(15) 



Following Abrikosov, we choose a general solution of 
Eq. ( 12 ) in a form of a sum over lowest Landau levels |21J : 



where 

hn{y,z,z) = exp|-^| 



Lb — \ 



(16) 



},(17) 



(18) 



is the magnetic length and v is an arbitrary parameter. 

We anticipate that the solution corresponds to a peri- 
odic lattice structure, and thus we assume that the vari- 
able n takes integer values in Eq. ([l2|, n G Z. The so- 
lution is parametrized by (arbitrary, in general) complex 
parameters Cn- In order to ensure a regular structure 
of the lattice, the coefficients Cn are usually chosen in a 
periodic manner: 



(19) 



The solution with N = 1 defines the square lattice (with 
all Cn's being equal). It corresponds to the original 
Abrikosov's solution [22 . 

In order to determine the optimal minimal-energy val- 
ues of the vortex lattice parameters (A/", Cn and v) in 
the GL model one usually solves classical equations of 
motion. In this article we would like to use an explicit 
energy minimization over these parameters because this 
approach is more suitable for a nonlocal energy functional 
of the p-meson sector of the QCD vacuum. 

The values of the coefficients Cn in Eq. ([l6| and the 
parameter v in Eq. (17) are determined by a minimiza- 
tion procedure of the action of the theory ([2| in strong 



magnetic field. Since the action density is a function of 
the transversal coordinates xi and X2, it it convenient to 
consider the energy density which is averaged over the 
transversal plane: 



-Bl 



(m2-eBext)(|</'|') 



(20) 



Using the explicit form of the solution given by 
Eqs. (16) and (17), one can show that 



M/2 

ni=-M/2 n2GZn3GZ 
^ni+n2^ni^ni+n3^ni+n2+n:i ' 



1 



1 



'I'^ = 2Mm-?ooM 



(21) 

2 / 2 I 2 \ 

(22) 



In the case of the A/'-fold symmetry (|19|) one gets: 
1 



]') = 

|4\ _ 



N-1 

E 

N-1 



Cn 



1 

2NM 

I I ni=0n2GZn3GZ 
^ni-\-n2^ni^ni-\-n3^ni-\-n2-\-n3 ' 



EEE 



where the iV-fold symmetry (19) is assumed. 



In the Bogomolny hmit of the couplings, 



A 



(23) 



(24) 



(25) 



the energy ( 20 ) depends only on the spatial average of the 
condensate squared At 5ext < the minimum 

of the energy density 



cmm 



1 



(25ext — Be) Be , 



(26) 



corresponds to the following value of the condensate: 

1 



2\ min 



-{Be - ^ext) • 



(27) 



Notice that in the Bogomolny limit an infinite amount 
of the lattice structures correspond to the same mini- 
mum of the energy functional (26) because the energy 



minimum is determined only by the specific value of the 
condensate (27). The mean condensate defines - accord- 
ing to Eq. (21 ) - a mean value of the squared coefficients 
Cn, and does not determine their precise values. There- 
fore, in the Bogomolny limit the lattice is no more regu- 
lar. The later fact nicely fits the observation that paral- 
lel Abrikosov vortices do not interact in the Bogomolny 
limit, so that a periodic lattice structure is unlikely to be 
formed. In the type-II regime ([5| the vortices repel each 
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According to Eq. (13), the magnetic field is concen- 



FIG. 1. Minimal-energy vortex lattices for N — 1 (left) and 
N = 2 (right) lattice symmetries in the GL model. 



trated at the positions of the vortices, so that the distri- 
bution of the magnetic field is visually similar to Fig. [l] 



III. CONDENSATION OF p MESONS IN QCD 

The conventional superconductivity is driven by the 
condensate of the Cooper pairs. The superconductivity 
of vacuum in a strong magnetic field is guaranteed by 
quark- ant iquark condensates which carry quantum num- 
bers of (charged) p mesons [14 . Below we consider the 
electrodynamics of p mesons in strong magnetic field. 



other so that the regular (periodic) vortex lattice is an 
anticipated structure of the ground state in this case. 

Since we are working with a type-II superconduc- 
tor ([5|, the last term in the energy density is always 
positive ( pQ| ). Therefore the global minimum of the en- 
ergy ( [2Q| ) - as a function of the parameter N which de- 
termines the symmetry in Eq. ([l9| - corresponds to a 
global minimum of a dimensionless quantity. 



,12^2 



which is known as the Abrikosov ratio |2T] 



(28) 



The 



Abrikosov ratio (28) in the ground state is independent 



of the value of the external magnetic field. 

For the simplest lattice of the square type. A/" = 1, 
the Abrikosov ratio ([28]) is f3{N = 1) = 1.180 which is 



reached at = 1. However, at A^ = 2 the Abrikosov ratio 
(and, consequently, the energy) reaches its global mini- 
mum, I3{N = 2) ^ 1.1596, at the following parameters: 



N = 2: Ci = ±iCo , 



4/o 

u = ^ ^ 0.9306 . (29) 
V2 



This minimal-energy periodic pattern corresponds to the 
equilateral triangular lattice (which is sometimes called 
"hexagonal" lattice). At higher odd values of the 
Abrikosov ratio is higher than the N = 2 minimum (for 
example, Pn=3 = 1.167) while at even values of N the 
minimization converges to the two-fold pattern corre- 
sponding to the triangular lattice {/3N=2k = Pn=2 with 
k £ Z). A nice review of the GL theory of type-II super- 
conductors in magnetic field, and, in particular, a review 
of the vortex lattice structure, can be found in Ref. [23] . 

The square {N = 1) and triangular {N = 2) Abrikosov 
vortex lattices are visualized in Fig. [l] The darker point- 
like regions indicate the positions of the cores of the vor- 
tices characterized by a reduced value of the condensate 
(/) = (j){xi,X2). The condensate vanishes at the center 
of each vortex, = 0, while the phase of the scalar 
field winds around a geometrical center of each vortex, 
arg^ = where is a two-dimensional azimuthal an- 
gle corresponding to a coordinate system centered at the 
vortex's origin. 



A. General properties 

1. Electrodynamics of p mesons 

The p meson is a vector (spin-1) particle made of a 
light (up or down) quark and a light antiquark. A self- 
consistent quantum electrodynamics for the charged and 
neutral p mesons is described by Djukanovic-Schindler- 
Gegelia-Scherer (DSGS) Lagrangian |24j: 



1 



1 



t2 



(30) 



.1 „(0) (0)m<^ , !!!p JO) {0)f. , 



which extends the vector meson dominance model [25^ 
with the Maxwellian U{1) sector by adding all allowed 
interactions of both charged, p^= p~ = {p\^^ — ip^^^)l\f2 

and p'^ = pI^^ and neutral, p^^\ mesons with the electro- 
magnetic field Ajj^. The tensor quantities in (30) corre- 
spond to various strength tensors. 



- d o^^ 



f(0) 

J LLV 



,(0) 



- d o^^^ 



,(0) _ fiO) 



Wsiplpu 



P^v 

PflU D^Py - DyP^ , 

where the covariant derivative is 
D 



Pfipt) ' 



(31) 
(32) 

(33) 
(34) 

(35) 



/i 

and gs = gpnTr ~ 5.88 is the pirir vertex coupling. 

The Lagrangian ( [3Q] ) respects the local U{1) symmetry: 

r pp{x) ^ e*-(-)p^(x). 



(36) 



2. The p-meson condensation due to strong magnetic field 



The last term of the DSGS Lagrangian (30) describes 



a nonminimal coupling of the p mesons to the electro- 
magnetic field which implies the anomalously large gy- 
romagnetic ratio {g = 2) of the charged mesons. In 
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the background of a strong enough magnetic field a spin- 
one particle with gyromagnetic ratio g = 2 should ex- 
perience a tachyonic instability towards formation of a 
Bose-Einstein condensate. 

In a simple way the condensation can be explained as 
follows [M]. A free charged relativistic spin-1 particle 
with the gyromagnetic ratio g = 2 and the mass m has 
the following energy spectrum in a background of an ex- 
ternal magnetic field ^ext = (O^O^^ext)' 

4,.. (Pz) = + (2n - 2s, + l)e5ext + • (37) 

where 5^ = — 1, 0, +1 is the spin projection on the field's 
axis i = ^3, n ^ is a nonnegative integer number 
(which, together with 5^, labels the Landau levels), and 
Pz is the particle momentum along the field's axis. The 
ground state (with n = 0^ = +1, =0) has the 
following (squared) energy: 



^o(^ext) 



m 



eB, 



(38) 



The lowest energy of the charged p-meson in the exter- 
nal magnetic field becomes purely imaginary quantity if 
the magnetic field exceeds the critical value Be = vn? je. 
Exactly the same result, Eq. (38), can be derived from 



the DSGS model (30) in a strong magnetic field back- 



ground [14 . The condensed state corresponds to the fol- 
lowing p- meson wavefunction: 



Pi 



-ip2 = p , Po = P3 = . 



Thus, at ^ext > the system experiences a tachyonic 
instability towards the Bose-Einstein condensation of the 
spin-1 particles. A similar condensation effect - a con- 
densation of the ly-bosons - is suggested to happen in 
the standard model of electroweak interactions in a back- 
ground with a much stronger magnetic field [26? ]. An- 
other example is represented by a pure Yang-Mills theory, 
in which the gluons may condense in a c/iromomagnetic 
field background [27 . 

The condensation of the p mesons in the vacuum takes 
place at the critical strength of the magnetic field ([T]) 
defined by the p- meson mass rrip = 775.5 MeV. The con- 
densation of the charged excitation - such as the p meson 
- should lead to an electromagnetic superconductivity of 
the vacuum. We would like to stress that the super- 
conducting state possesses the quantum numbers of the 
vacuum (i.e., the local electric charge density is zero at 
every point of the space-time). 



3. Energy density of the ground state 

Let us consider the spontaneous condensation of the p 
mesons in the background of the magnetic field (|4| in the 
vicinity of the phase transition: B > Be with B—Bc<^Bc. 
Since the external field is taken to be slightly stronger 
than the critical value ([T]), the p- meson condensate is 
small, IpI <C m^, so that the equations of motion for p- 
meson electrodynamics (|3Q|) can be linearized. Following 



the example of the Ginzburg-Landau model described in 
Section |n| we consider static xs-independent solutions 
which may, however, be inhomogeneous in the transversal 
(xi, X2) plane. 

The classical equations of motion in the magnetic field 
background are quite complicated and we refer an inter- 
ested reader to Ref. 14 for a detailed discussion. It turns 



out that the equation for the p- meson condensate ( 39 ) in 
the overcritical magnetic field {B > Be) is similar to 
the one for the Cooper pair condensate in the subcritical 
magnetic field {B < Be) in the Ginzburg-Landau model: 



T)p={d--B,^tz)p = 0, 



(40) 



where the covariant derivative is given in Eq. ([nj and 
the complex wavefunction p = p{z) is defined in Eq. (39) 
Solutions of Eq. (40) are inhomogeneous in the 



transversal (xi,X2) plane, and these inhomogeneities in- 
duce condensation of the neutral mesons, p'^^^ 



• (0) 

^P2 ' 



with pQ 



(0) 

Pi 



P^'Hx^) 



2igs 



-di 



rd\p\' 



(0) _ ^(0) 



0. In Eq. (41) dl = d"^ 



dl 



-dl 



— 2ix±) = —Ko{m\x±\) , 



(41) 



(42) 



(39) is the two-dimensional Euclidean propagator of a scalar 



particle with the mass of the neutral p^^^ meson, 



mo 



m, 



(43) 



and Kq is a modified Bessel function. 

The magnetic field is also an inhomogeneous function 
of the transversal coordinates. 



B{x^) = B, 



2emQ 



-dl 



mf, 



\p{z)\'-{\p\'- 



(44) 



where the last term, defined in Eq. (14), guarantees the 
conservation of the magnetic flux (15). Equation (44) is 
a p- meson analogue of the GL relation ([T3|. 

Notice that in the vacuum subjected to the strong mag- 
netic field, the neutral condensate (|4T|) and the magnetic 



field (44) depend on the (charged) p-meson condensate 
nonlocally contrary to the local relation between the mag- 
netic field and the Cooper pair condensate ( 13 ) in the GL 
model for the ordinary superconductivity. 

The energy density of the vacuum in the presence of 
the p-meson condensate is given by the following formula: 



(Too) = A + 2K - eSext)(|pn + 2e'{\p\') 



2\2 



were T^jj is the energy- momentum tensor. 



(45) 



(46) 
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corresponding to the DSGS model (30), and the brackets 



(. . . ) indicate the average over the transversal (xi,X2) 
plane (14). Notice that contrary to the GL model (20), 



the energy density of the p- meson condensate (45) con- 



tains a nonlocal positive definite {gg ^ e) quartic term. 

B. Physical properties of the ground state 

1. Superconductor vortices: the lattice structure 

Following our experience in the GL model we repre- 
sent the solution of the p-meson condensate in a manner 
similar to Eq. ([l6|, 

p{z) = J2Cr.hn{l^,^,^), (47) 



where the function hn is given in Eq. (17) and Lb is the 



magnetic length ( |18| ). The coefficients are assumed to 
obey the A/'-fold symmetry (19). 

The quadratic term in the energy density (45) 



can be evaluated with the help of E g. (23[ ), while the local 
quartic term is absent in Eq. (p5). Instead, the con- 
densate p is stabilized by a nonlocal quartic term, which 
is proportional to the following nonlocal functional: 



Q[p] = ( \p\ 



-A + ml 



\P\' 



(48) 



which can conveniently be represented as follows: 
1 f (Pk w? 



Q[p\ 



Area_L J {2ttY + 



\q{k;p)\\ (49) 



where 



q{k;p)= / dxi / dx2e^^^"^+^^^"^|p(xi,X2)|', (50) 



with p) = q{-k; p) and X2) = p{xi + 1x2)^ 

Let us shift the momentum, k k/Ls^ and the coor- 
dinate, z Lbz, so that the new variables k and z are 
now dimensionless quantities. Then quantity (50) takes 
the following form 

q{k;p) = LlY^ ^ C„,C:^F„„„,(fci,fc2), (51) 

niGZ n2GZ 

where 

Hn,,n2{ki,k2) = J d^ze'^-^ hn,{z,z)hn2{z,z) 

= V2ii8(k2 - 2iiv(ni - n2))e-^^'(^^-^2)' 
ni + 712 k\ 



exp < ikiv- 



with k • z = kiXi + k2X2- 



8^1' 



(52) 



Next, we substitute Eq. (51) into Eq. (49), 



Area^ 
1 



S — 1 1 



) -f^ni,n2 (^)-^n3,n4 ( ^) ? (^3) 



(27r)2 k^ + m? 
and take into account the following relation: 

5(^2 — 27ru{ni — 712)^ s(^k2 + 27Tu{ns — 114^ 



)) 



N2 



■n2+n3 — n4 



,0 5 [k2 - 27iv{ni - 712)) , (54) 



where N2 is the number of elementary cells in X2 direction 
of the Abrikosov lattice [notice that 5{x) is the Dirac 5 
function while ^^,0 is the Kronecker symbol]. Next, we 
make the following shifts of the integer- valued variables: 
Til ^ Til + 7i2, ^3 ^ na + 77,4, and 77,4 ^ 77,1 — 712 + ^3. 
The result is 



Q[p] = 



1 N2L% 

2 •'"^."3 



EEE/» 

niGZn2GZn3GZ 



^ni+n2^ni'-^ni+n3^ni+n2+n3 ' 



(55) 



where the matrix fn2,n3 is given by the following relation: 

/n2,n3(/^5 ^) = rr / 7.2 , /i^2,.2^2 , ..2 ' (^^) 



27r k^ 



p- 



We expect that in the infinite-mass limit, 777,0^00, 
the functional (48) should converge to the local quantity. 



hm Q[p] 

mo^oo 



And, indeed, in the limit p = ttiLb 
the matrix (56) the following expression: 



(57) 

00 one gets for 



lim /n2,n3(M) 



2 / 2 I 2 \ 

-TTU (n2+n3) 



/W\ 



(58) 



SO that the nonlocal functional (59) converges to the local 
expression (24), as expected. The function /n2,n3 is not 



symmetric under discrete 7r/2 rotations in the (77,2,723) 
plane at small values of the dimensionless parameter p. 
However, as the mass parameter p increases, the function 
/ approaches the symmetric Gaussian function (58). 



Finally, taking into account the A/'-fold periodicity ( 19) 



of the coefficients C, one can rewrite Eq. (55) as follows: 

I AT-l 

Q[p] 



2Ar ^ ^ 

ni=0 n2GZ 
^ni-\-n2^ni^ni-\-n3^ni-\-n2-\-n3 ' 



2N ^ ^ ^ /n2,n3(^0^B,i^) 

ni=0 n2GZ n3GZ 



(59) 



The structure of the energy functional (45) indicates 
that an analogue of the Abrikosov ratio (28) in the case 
of the p-meson condensation is as follows: 



/ \P? 



\P\' 



\(H2)-A + m2(H2) 



Q[p] 

/U|2\2 



(60) 



N 




MeV 


V 


Co, MeV 


Ci 


C2 


Cs 


C4 


type 


1 


1.02393 


26.7334 


1 


10.9516 










square 


2 


1.01920 


26.7656 




10.5977 










hexagonal 


3 


1.02056 


26.7567 


0.95581 


10.9744 


Co e'^^ 


Co e'"^^ 






parallelogrammic 


4 


1.01920 


26.7656 


^/V2 


10.5977 


Co 


-Co 


-Co 




hexagonal 


5 


1.01954 


26.7633 


0.93883 


10.9773 


Co e-^^'^^ 


Co e-^^^^^ 


Co e-^^'^^ 


Co 


parallelogrammic 



TABLE I. The parameters of the p-meson lattices at Bext = l.OlBc: the dimensionless ratio /3p, Eq. (|60|, the condensation 
energy per unit volume, Eq. (61), and the structure constants Cn for the first five lattice types N — 1,...,5 (here 
— 2n/k). Examples of certain lattices are shown in Fig.|2] The global minimum in energy is reached N — 2. 
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FIG. 2. The density plots of the superconducting density in the p-meson vortex lattices with minimal energy at fixed lattice 
type (a) iV = 1, (b) iV = 2, (c) = 3, and (d) = 5 at B = l.OlBc. The corresponding magnetic length (flSl is Lb = 0.634 fm. 



The darker regions corresponds to positions of the p vortices where the superconducting density is suppressed. 



The minimum of the energy functional corresponds to the 
minimum of the new dimensionless parameter Pp. Con- 



trary to the Abrikosov ratio (28), the quantity (60) de- 
pends on the strength of the magnetic field B. In the "lo- 
cal" (and, unphysical) limit mo oo the quantity (60) 



is reduced to the Abrikosov ratio ( 28 ) . 

We minimize numerically the mean energy density ( 45 ) 
as a function of (generally, complex) lattice parameters 
Cn, n = 1, . . . , A" for a fixed value of A" = 1, . . . , 8. We 
have found that the condensation energy. 



5£ = {£) 



1 



(61) 




1.00 1.05 1.10 1.15 1.20 1.25 1.30 

B/B, 



reaches its minimum at the equilateral triangular lattice 
with N = 2 and Cq = iCi, similarly to the case of the 
Abrikosov lattice in the GL model. All lattices with odd 
values of A" possess higher energies while all even-A^ lat- 
tices are reduced to the N = 2 case. 

In Table [l| we present the results of a numerical min- 
imization of the energy functional (45) for fixed values 
of A" at B = 1.015c. As in the case of the GL model 
for a type-II superconductor, the difference in minimal 
condensation energy 6£ for the different values of A^ is 
very small. The global minimum in SS and in the pa- 
rameter Pp is reached for even values of A" which corre- 
sponds to the same equilateral triangular lattice which 
is sometimes called "hexagonal lattice" . Notice that the 
explicit parameterization of the parameters Cn at differ- 
ent (even) values of A" may not correspond to multiple 



FIG. 3. The parameter /3p, Eq. (60), for different types of 
lattices. A" = 1, 2, 3, as a function of the magnetic field B. 



repeated copies of the N = 2 solution (this fact is clearly 
seen for the case of A^ = 4). The simplest square lattice is 
realized at A" = 1 while odd values of A" ^ 3 correspond 
to a general (called "oblique" ) type of the 2 dimensional 
lattice. 

Examples of the minimal-energy lattices at A^ = 
1,2,3,5 formed in the magnetic field B = 1.01 Be are 
presented in Fig. [2] The minimal-energy lattices with 
A" = 4, 6, . . . coincide with the N = 2 solution. 

In Fig. Is] we show the parameter Eq. (60), as a 
. oft 



function of the magnetic field B in the vicinity of the 
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critical magnetic field Be. The minimum of the quantity 
f3p - and, as a consequence, the minimum of the energy 
density (45) - are reached at A/" = 2 for all studied val- 



ues of the magnetic field. Due to specific values of the 
phenomenological parameters of the DSGS model (30) 



- which describes the electrodynamics of the p-meson 
excitations in the QCD vacuum - the difference in en- 
ergies between visually different lattices (c./. Fig. [2| is 
tiny. For example, at B = 1.01 5c the difference in the 
condensation energies between the square, TV = 1, lat- 
tice and the equilateral triangular, = 2, lattice is less 
than 0.5%. The relative difference in the corresponding 
dimensionless P parameters is of the same order. 

In the first study of the vacuum superconductiv- 
ity in Ref. |^ the analysis of the ground state was 
done in assumption - following the pioneering work of 
Abrikosov [22] - that the vortex lattice has a square pat- 
tern with N = 1. However, the more detailed analy- 
sis of this paper indicates that the real ground state of 
the superconducting vacuum is a triangular lattic^with 
N = 2. Since the difference in most important bulk pa- 
rameters (e.g., average energy, mean conductivity etc) 
between the square lattice and its possible conformations 
is very small, the square lattice is a very good approxi- 
mation for calculation of the bulk properties of the real 
vacuum state. In this article we explore the correct (tri- 
angular vortex) state of the vacuum superconductivity 
which is important for local properties of this unusual 
phase. 



2. Superconducting condensate and energy density 

From now on we concentrate on the triangular lattice 
which corresponds to two-fold symmetric {N = 2) pa- 
rameter spac e (p^O] ), with the lowest-energy parameters 
given by Eq. (|29[). 

In Fig . El (top) we show the mean condensation energy 
density ( [M] ) as a function of the magnetic field. One can 
clearly see that at B < Be the condensation energy is 
zero while it becomes negative dX B > Be due to conden- 
sation of the charged p mesons. In order to characterize 
the later property, we notice that for the equilateral tri- 
angular vortex lattice the mean squared superconducting 
condensate is related to the coefficient Co of the solu- 



tion (47) as follows: 



(fx\p{x)[' 



2^ 



Col 



(62) 



In Fig. |4] (bottom) we plot the square root of the mean 
of the squared condensate ( 62 ) as a function of the mag- 



netic field. It is clear that Sit B > Be the superconducting 
state with a nontrivial condensate p ^ is energetically 




^ Soon after the Abrikosov's paper [22| the vortex pattern in a 
type-II superconductor was shown to be triangular [21| . 



FIG. 4. At B > Be the superconducting state is more energet- 
ically favorable compared to the trivial vacuum state: (top) 
the condensation energy (61) becomes negative due to emer- 
gence of the superconducting condensate (|62|, \p\ = 
(bottom) at the critical magnetic field B — Be with Be given 
in Eq. ([T]). The quantities are shown for the minimal-energy 
(equilateral triangular) lattice structure. 



more favorable compared to the trivial vacuum state with 
p = 0. The stronger magnetic field the larger the gain in 
energy due to the condensation of the p mesons. 

The vortex structure of the superconducting ground 
state can easily be seen given the behavior of the su- 
perconducting order parameter p. The p condensate has 
a characteristic form p{z^ z) oc z — zq near the position 
^0 = ^1,0 + ^^2,0 of an elementary vortex. Thus, at the 
vortex core z = zq the condensate should vanish as a lin- 
ear function. Moreover, in the local vicinity of the point 
Zq the phase cpp = arg p of the condensate should behave 
as a polar angle. Thus, the phase (pp winds around the 
position of the vortex and the winding corresponds to a 
topological stability of the vortex. Since the phase is de- 
fined modulo 27r, the phase (fp should experience cuts in 
the (xi,X2) plane. At branches of these cuts the phase 
(Pp experiences the quantized jumps (pp ^ (pp± 27r. 

The absolute value and the phase of the superconduct- 
ing condensate are shown in Fig. [5] and Fig. [6j respec- 
tively. From Fig. |5] one can see that the condensate van- 
ishes linearly {p oc \z — zo\) around a set of isolated points 
which form the equilateral triangular pattern. Fig. |2|b). 
The phase of the condensate. Fig. [6j experiences jumps 
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Ipl, 10 
MeV 5 




FIG. 5. Absolute value of the superconducting condensate p, 
Eq. (471 at B = 1.01 Be in the transversal (2:1,2:2) plane. 




FIG. 6. The density plot o f th e phase cpp = argp of the 
charged yo-meson condensate ( 47 ) as a function of transversal 
coordinates xi and X2. The white lines corresponds to the 
cuts in the phase, and the variations in color illustrate the 
behavior of the magnitude of the phases. The endpoints of 
the cuts mark the positions of the superconductor (p) vortices. 
A three-fold difference in coordinate scales of this figure and 
Fig. [5] should be stressed. 



at one-dimensional manifolds (semi-infinite lines) which 
start at the points where the condensate p is vanishing, 
and end at spatial infinity. The position and the shape 
of the cuts can be changed by the U{1) gauge transfor- 
mations (36), while the endpoints of the cuts are gauge 



invariant quantities. 

Thus, the p = points of the charged condensate do in- 
deed mark the positions of the superconductor vortice^ 
and organize themselves into the equilateral triangle lat- 
tice. 



3. Electric currents and superconductivity 

The electric current density of the p meson degrees of 
freedom can be derived from Eq. (ISOl): 



^/i = ^e[p^Vr^/i - p'^pU + d^'iplp^ - pIp^)] 



(63) 



In the ground state the longitudinal components of the 
electric current are vanishing, Jq = J3 = 0, while the 
transversal current, 



J^{x^) = Ji{x±) + iJ2{x±) , 



(64) 



is a nonlocal function of the superconducting conden- 
sate [H]: 

J±{x±) = 2ieml ■ ( _g2 ^^2 I/^P)(^±) . (65) 

The nonlocal nature of the relation between the transver- 
sal electric current (65) and the charged condensate p 



distinguishes the QCD vacuum from the GL supercon- 
ductor ([2|. 



The electric current ( 65 ) is a persistent current of the 



charged p-meson degrees of freedom. The current origi- 
nates from the quarks and antiquarks which popup from 
the virtual state and form a condensate with the quan- 
tum numbers of the charged p-meson. This current is 
always present in the superconducting phase B > Be 
and it vanishes in the normal phase of the vacuum. 

In order to simply further numerical calculations it is 
useful to represent the basic nonlocal part of Eq. ( 65 ) 



-di 



(66) 



as follows: 



R{x±;p) = J 



cPk 



q{-k;p),{67) 



where the quantity q{k] p) is defined in Eq. (50) and p 
is given in Eq. (47) with N = 2 symmetry. For the 



triangular lattice with Ci = iCo one can explicitly show 
that 



^ We call the topological defects in the charged p condensate as ^^^-^ (0) condensate 

"superconductor vortices" in order to distinguish them from "su- 
perfluid vortices" which are similar vortexlike defects in the neu- 
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niGZn2GZ\ 



SO that the nonlocal quantity (66) is 

\Co\' 



2W\ 



ni GZ n2 G2 



TT 

"21 



-2 2 I 2 2 



[i-(-i)"-](-i)"^"n-i)(i-"-)/V _^.,^/; 

47r2[i/-2(ni + l/2)2 + jy2n2] +^2 e 



(68) 



i^/ci — 27r^77,i + J(/c2 — 27rz^77,2) exp| — ^ ^(^^^2) ^^^^2 I 



(69) 
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FIG. 7. The electric currents of the charged quark con- 
densates in the ground state lattice at B = 1.01 Be in the 
transversal (xi,X2) plane. The reddish (grayish) areas corre- 
sponds to the stronger (weaker) current. 



with ji = tyi^Lb^ where the magnetic length is given 
in Eq. ([T8| . This expression is very convenient for numer- 
ical calculations because the sums in Eq. (69) converge 



very quickly. The derivatives of Eq. ( 69 ) with respect to 
the transverse coordinates xi and X2 can be taken quite 
straightforwardly. 

In Fig. [7] we show the transverse electric current ( 65 ) 
in the ground state of the superconducting vacuum at 
the magnetic field B = 1.01 All current lines are 
organized in the hexagonal pattern. The center of each 
hexagon corresponds to the position of a superconductor 
vortex. 

In transversal directions the strong magnetic field con- 
fines the local electric currents of charged condensates to 
the hexagonal cells. The size of cells are of the order of 
the lowest Landau level of the p-meson condensate, that 




x-i , f m 



FIG. 8. The superconductivity coefficient Eq. (72), in the 



modified London law (70), is shown as a function of xi and 
X2 transverse coordinates at the magnetic field B = 1.01 Be. 



is of the order of the magnetic length ( |18| ) {Lb ^ 0.6 fm 
dit B ^ Be). A relatively week (|^| <C \B\) electric field 
E ± B cannot create a transversal ( "intra-cell" ) electric 
current because such a current would involve an excita- 
tion at a first Landau level which is separated from a 
lowest Landau l evel by a large energy gap of the order 
of El — Eq ^ ^/\eB\. Thus, the global electric currents 
are drastically suppressed in the transverse directions. 
Therefore, in the transverse direction the vacuum state 
behaves as an insulator. This is the very reason why the 
Meissner effect is absent in the superconducting ground 
state [14J so that the emerging superconductivity does 
not screen the external magnetic field. 

However, the longitudinal electric currents can exist 
as the motion of the quarks along the axis of the mag- 
netic field is not restricted. If one applies a weak electric 
field E = (0,0, £^3) parallel to the strong magnetic field 
B = (0, 0, B) then the electric currents - induced by the 



11 



external field - satisfy a modified London equation [14] : 

SJ3(x) 9Jo(x) 



(70) 



dxo dxs 

while the induced transverse electric currents are obvi- 
ously zero, so that the transverse current J± is unmodi- 
fied by the external electric field directed along the mag- 
netic field axis, Eqs. and (|65|, and 



dJk{x) dJ^{x) 



dx. 



dxk 



0, 



(71) 



with = 0, . . . , 3 and A: = 1, 2. 

The London-type equation is typical for superconduct- 
ing systems as it characterizes a conducting state without 
resistance. The superconductivity is characterized by the 
quantity n = hz{x±)^ which is a nonlocal function of the 
superconducting condensate: 

n{x^)=4e^ml-(—J—^\p\^){x^). (72) 

In the superconducting ground state the "supercon- 
ducting transport coefficient" has a hexagonal lattice 
structure. Fig. [S] Due to nonlocal nature of the rela- 
tion between the transport coefficient and the super- 
conducting condensate (72), the superconductivity is not 



completely suppressed inside the vortices. This property 
distinguishes the vacuum superconductivity from an or- 
dinary one. In the GL model of superconductivity, the 
mentioned relation is local and the suppression of the 
condensation in the center of the Abrikosov vortex im- 
plies strong suppression of the superconductivity by the 
vortices. 

Thus, at high magnetic field the vacuum becomes an 
anisotropic inhomogeneous superconductor. The super- 
conductivity is anisotropic since the vacuum can super- 
conduct along the direction of the magnetic field only 
and in the transversal directions the vacuum behaves as 
an insulator. The superconductivity is inhomogeneous 
due to the coordinate dependence of the superconduct- 
ing transport coefficient k,^ Eq. (72), Fig. [s] 



4. Neutral condensate and superfluid vortices 

It is known that the condensation of the charged (su- 
perconducting) field pjj^ leads to the induced condensation 

of a neutral, superfluid-like field p^^^ in the ground state 
of the vacuum [14 . The longitudinal components of the 

neutral condensate are zero, p'q^ = p^^ = 0, while the 
transverse components of the neutral condensate, are, in 
general, nonvanishing: 



P^°\x^) = 2ig, ■ ( _g2+^2 H')(^^) ' (73) 



where the quantity p*^^^ = pl ^ -\- ip2 is a complex field. 

The superfluid current - which can be interpreted 
as a flow of the neutral condensate - can be derived in a 




Xj , f m 



FIG. 9. Absolute value of the superfluid condensate p''^\ 
Eq. (73), at B = 1.01 i^c in the transversal (xi,X2) plane. 



similar manner to the electromagnetic current (63) via 



variation of the action of the p-meson electrodynam- 
ics (30) with respect to the neutral field p^^\ In the 
ground state of the vacuum the longitudinal superfluid 
current vanishes, Kq = Ks = 0, while the transversal 
current. 



K±{x±) = Ki{x±) -^iK2{x±) , 



(74) 



is a nonlocal function of the superconducting condensate: 

In Fig.|9]we show the absolute value of the neutral con- 
densates p^^^ as a function of the transverse coordinates 
xi and X2 for a slightly overcritical background magnetic 
field, B = 1.01 One can compare this figure with 
Fig. |5] where the charged condensate is plotted for the 
same set of parameters. Firstly, we notice that the neu- 
tral condensate is much weaker compared to its charged 
counterpart, Ip*^^^! <C |p|. Secondly, the geometrical pat- 
tern of the neutral condensate is much more involved 
compared to the charged one. In particular one observes 
that the neutral p-meson condensate vanishes in a denser 
set of points compared to its superconducting counter- 
part. The latter set of points is organized in a triangular 
pattern. Fig. |2j As the phases of the charged p field are 
winding around these points (Fig. |6|, we concluded that 
these points correspond to the superconductor vortices 
("p vortices"). Following this analogy, one can suspect 
that the points of vanishing neutral condensate may cor- 
respond to "superfluid" vortices and antivortices. In or- 
der to figure out the topological structure of the neutral 
condensate we should check whether the phase of the 
neutral p^^^ field winds around these points or not. 



In Fig. [To] we show the density plot of the phase of 
the neutral p-meson field (73) in the transverse plane 
(xi,X2). The white lines corresponds to the cuts in the 
phase, so that the phase of the neutral condensate winds 
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FIG. 10. The density plot of the phase of the neutral p-meson 



field (73) as a function of transverse coordinates xi and X2- 



The white lines corresponds to the cuts in the phases, and the 
variations in color illustrate the behavior of the magnitude of 
the phase. 



around the endpoints of these lines. At these endpoints 
the absolute value of the neutral field vanishes and the 
phase becomes undefined. Thus, the endpoints corre- 
spond to the (superfluid) vortices in the neutral p-meson 
field. Surprisingly, the superfluid vortices and antivor- 
tices comes always in pairs so that the net vorticity of 
the neutral field is zero. 



The superfluid current (75) circulates around the su- 



perfluid vortices and antivortices in clockwise and coun- 
terclockwise in a manner which is visually similar to 

Fig. [3 

Finally, we would like to mention that the transversal 
electric current = Ji + zJ2, Eq. (65), the supercon- 
ducting transport coefficient Eq. (72), the neutral con- 
densate p^^^ = p^^^ + ^P2^\ Eq. (73), and the superfluid 



current K_\_^ Eq. (75), can be mutually related to each 
other: 



emf{ 



Ws 



(76) 
(77) 



where d = di -\- 182 and d\ = d1 ^ 82 • This relation 
is valid in the leading order in the vicinity of the phase 
transition, B > Be. 



0.0 




a superfluid vortex on top 
of a superconductor vortex 



superfluid antivortex 



FIG. 11. The periodic vortex structure of the vacuum ground 
state is superimposed on the density plot (shown in shades 
of the green color) of the absolute value of the neutral p- 
meson condensate ( [73| at the magnetic field B — 1.01 Be- 
Each superconductor vortex [the large red circles] is always 
superimposed on a superfluid vortex (the small blue disks 
marked by the plus signs) forming an equilateral triangular 
lattice. The isolated superfluid vortices and antivortices (the 
small yellow disks with the minus sign) arrange themselves in 
the hexagonal lattice pattern. 



on the density plot of the absolute value of the neutral 
meson field (73). The vortex pattern is quite remarkable: 



• The superconductor vortices organize themselves in 
an equilateral triangular lattice; 

• The superfluid vortices and antivortices are orga- 
nized in a honeycomb (hexagonal tiling) pattern; 

• A center of each hexagon is occupied by a super- 
conductor vortex; 

• Each superconductor vortex is always superim- 
posed on a superfluid vortex; 

• One superconductor vortex is accompanied by 
three superfluid vortices and three superfluid an- 
tivortices so that the net superfluid vorticity is zero. 



IV. CONCLUSION 



5. Superconductor and superfluid vortex lattices 

We show both the superconducting and superfluid vor- 
tices in Fig. 11 The vortex locations are superimposed 



In this paper we determine the structure of inhomo- 
geneities of the quantum vacuum in a superconducting 
state. The superconducting state of the vacuum is real- 
ized at a strong magnetic field background if the strength 



13 



of the field exceeds the critical value ([T]), Ref. [14]. The 
superconductivity is associated with a spontaneous emer- 
gence of quark- ant iquark condensates which carry quan- 
tum numbers of charged p mesons. The p-meson con- 
densate in the ground state turns out to be an inhomo- 
geneous structure made of the vortexlike defects in the 
p-meson field. These superconductor vortices are paral- 
lel to the magnetic field axis, and in this paper we show 
that in the transversal plane they organize themselves in 
a equilateral triangular lattice. 

The condensation of the charged p mesons induces a 
(much weaker) superfluid-like condensate with quantum 
numbers of the neutral p^^^ mesons [14 . In this paper 
we have shown that the neutral condensate contains vor- 
texlike defects as well (by analogy we call them as the 



"superfiuid vortices"). We show that each of these su- 
perconductor vortices is accompanied by three superfiuid 
vortices and three superfiuid antivortices made of the 
neutral p meson condensate. Moreover, one of these su- 
perfiuid vortices overlaps with a superconductor vortex, 
so that each superconductor vortex is always accompa- 
nied by the superfiuid vortex. The superfiuid vortices 
organize themselves in a honeycomb lattice. 
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